On Darboux-Backlund Transformations for the 
Q-Deformed Korteweg-de Vries Hierarchy 



Ming-Hsien Tu^, Jiin-Chang Shaw^ and Chin-Rong Lee^ 
^ Department of Mathematics, National Chung Cheng University, 
Minhsiung, Chiayi 621, Taiwan 
^ Department of Applied Mathematics, National Chiao Tung University, 

Hsinchu 300, Taiwan 
^ Department of Physics, National Chung Cheng University, 
Minhsiung, Chiayi 621, Taiwan 
(February 9, 2008) 



Abstract 

We study Darboux-Backlund transformations (DBTs) for the g-deformed 
Korteweg-de Vries hierarchy by using the g-deformed pseudo differential op- 
erators. We identify the elementary DBTs which are triggered by the gauge 
operators constructed from the (adjoint) wave functions of the hierarchy. Iter- 
ating these elementary DBTs we obtain not only g-deformed Wronskian-type 
but also binary-type representations of the tau-function to the hierarchy. 

Mathematics Subject Classifications(1991): 58F07, 35Q53 
Key words: KdV hierarchy, q-deformations, Darboux-Backlund transforma- 
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Typeset using REVTgX 



I. INTRODUCTION 



There has been remarkable interest in the theory of sohton ever since the discovery 
of the Inverse Scattering Method (ISM) for the Korteweg-de Vires(KdV) equation(see, for 
example, fl]]). Although the ISM has been extended to many nonlinear systems which 
describe phenomena in many branches of science, the KdV equation still plays an important 
role in the development of modern soliton theory (see, for example 0]). In particular, many 
concepts were established first for the KdV equation and then generalized to other systems 
in different ways. 

A convenient approach to formulate the KdV hierarchy relies on the use of fractional- 
power pseudo-differential operators associated with the scalar Lax operator L = d"^ + u 
which was further generalized by Gelfand and Dickey |]^ to the A^-th KdV hierarchy which 
has Lax operator of the form = + UN-id^~^ -|- ■ ■ ■ + mo- In the past few years, there 
are several works concerning the extensions of the KdV hierarchy in Lax formulation, such 
as Drinfeld-Skolov theory and supersymmetric generalizations etc. The common 
features of these extensions show that they preserve the integrable structure of the KdV 
hierarchy and contain the KdV as reductions in some limiting cases . 

Recently, a new kind of extension called g-deformed KdV (g-KdV) hierarchy has been 



proposed and studied pHT^- In this extension the partial derivative d is replaced by the 



g-deformed differential operators (g-DO) dq such that 

(a,n^)) = lMjiI<fl (1.1) 

x{q — 1) 

which recovers the ordinary differentiation {dxf{x)) as q goes to 1. Even though many 
structures of the g-KdV hierarchy such as infinite conservation laws 0, bi-Hamiltonian 
structure 0, tau-function |TI|-|I1|, Vertex operators Virasoro and VT-algebras etc. 



have been studied however the g-version of Darboux-Backlund transformations (DBTs) for 
this system are still unexplored. It is well known that the DBT is an important property to 
characterize the integrability of the hierarchy Thus it is worthwhile to investigate the 
DBTs associated with the g-KdV hierarchy. Once this goal can be achieved, it will deepen 
our understanding on the soliton solutions of the hierarchy. 

Our paper is organized as follows: In Sec. II, we recall the basic facts concerning the q- 
deformed pseudodifferential operators (g-PDO) and define the A^th g-deformed Korteweg-de 
Vries (q-KdV) hierarchy. In Sec. Ill, we construct the Darboux-Backlund transformations 
(DBTs) for the A^th q-KdV hierarchy, which preserve the form of the Lax operator and the 
hierarchy flows. Iteration of these DBTs generates the g-analogue of soliton solutions of the 
hierarchy. In Sec. IV, the case for A^ = 2 (g-KdV hierarchy) is studied in detail to illustrate 
the g-deformed formulation. Concluding remarks are presented in Sec. V. 

II. Q-DEFORMED PSEUDODIFFERENTIAL OPERATORS 

Let us define the g-shift operator 9 as 

9{f{x)) = f{qx) (2.1) 
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then it is easy to show that 6 and dq do not commute but satisfy 

{d,e\f)) = q'e'idj), k e Z (2.2) 
Using (|1.1| ) and (|2rT| ) we have (dgfg) = {dqf)g + 9{f){dqg) which imphes that 

dqf={dqf)+e{f)dq (2.3) 

We also denote as the formal inverse of dq such that dqd~^f = d~^dqf = f and hence 



(2.4) 



fe>0 



In general one can justify the following g-deformed Leibnitz rule: 



dqf = E 



n 



ne Z 



(2.5) 



fe>0 L"'J q 

where we introduce the g-number and the g-binomial as follows 



[n\ 

n 
k 



-I g 



q-1 

[n]q[n -l]q---[n-k + l\^ 



(2.6) 
(2.7) 



For a g-PDO of the form 



i=—oo 



it is convenient to separate P into the differential and integral parts as follows 



(2.8) 



(2.9) 



t>0 



i<-l 



and denote (-P)o as the zero-th order term of P. The residue of P is defined by res(P) = n_i 
and the conjugate operation "*" for P is defined by P* = J2i{d*yui with 



d; = -dqO-' 

Then a straightforward calculation shows that 

(PQ)* = Q*P* 

where P and Q are any two g-PDOs. Finally, for a set of functions /i, /2, 
g-deformed Wronskian determinant Wq[fi, /2, ■ ■ ■ , /n] as 



(2.10) 

(2.11) 
/„, we define 



/i 

(dqh: 



fn 



(2.12) 
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With these definitions in hand, we have the following identities which will simplify the 
computations involving compositions of g-PDOs. 
Proposition 1: 

(p*)^ = {p^)\ {p*)_ = {p_y (2.13) 

id,'P)- = d;\no + d^'P- (2.14) 

iPd;')^ = iP)od;' + P-d^' (2.15) 

res(P*) = -^"^(res(P)), (5,res(P)) = res(9,P) - e(res(P(9g)) (2.16) 

res(P9-i) = (P)o res(9,-^P) = e-\{P*)o) (2.17) 

Tes{d;'PiP2d;') = Tes{d;\P;)oP2d;') + res(9-iPi(P2)o5;') (2.18) 

where Pi = (Pi)+ and P2 = (P2)+. 

The iV-th g-KdV hierarchy is defined, in Lax form, as 

9t„L=[P„,L], n = 1,2,3,- ■■ (2.19) 

with 

L = d^ + + ■ ■ ■ + no, P„ = L!^/^ (2.20) 

where the coefficients Ui are functions of the variables {x,ti,t2, ■ ■ ■) but do not depend on 

(tN, t2N, tsN, ■ ■ ■)• 

In fact, we can rewrite the hierarchy equations (|2.19| ) as follows: 

dt^B^n - dt^Bra ' [Pn, B^] = (2.21) 

which is called the zero-curvature condition and is equivalent to the whole set of equations 
of (|2.19|) . If we can find a set of functions {ui, i = 0, 1, ■ ■ ■ , A^ — 1} and hence a corresponding 
Lax operator L (or P„) satisfying ( ^.19] ) (or ( ^^.21] )), then we have a solution to the iV-th 
g-KdV hierarchy. 

For the Lax operator ( p.20|) , we can formally expand L^^^ in powers of dq as follows 

L^/^ = dg + so + sid-^ + ■ ■ • (2.22) 

such that (L^/^)^ = L which gives all the Si being g-deformed differential polynomials in 
{ui}. Especially, for the coefficient of d^~^ we have 

= So + e{s^) + ■■■ + e''-\so). (2.23) 

The Lax equations ( 2.19| ) can be viewed as the compatibility condition of the linear 

system 

L(j) = Ac/), dt„(j) = (P„0) (2.24) 

where (j) and A are called wave function and eigenvalue of the linear system, respectively. On 
the other hand, we can also introduce adjoint wave function ip which satisfies the adjoint 
linear system 

L*^ = /i^, ai> = -(P:^) (2.25) 

For convenience, throughout this paper, (pi {ipi) will stand for (adjoint) wave functions with 
eigenvalues Xi{fii), respectively without further mention. 
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III. ELEMENTARY DARBOUX-BACKLUND TRANSFORMATIONS 



In this section we would like to construct DBTs for the A^-th g-KdV hierarchy. To attain 
this purpose, let us consider the following transformation 

L L(^) = TLT-^ (3.1) 

where T is any reasonable g-PDO and denotes its inverse. In order to obtain the new 
solution {L^^^) from the old one (L) the gauge operator T can not be arbitrarily chosen. It 
should be constructed in such a way that the transformed Lax operator L^^^ preserves the 
form of L and satisfies the Lax equation (p.l9|) . From the zero-curvature condition (|2.21|) 
point of view, the operator should be transformed according to 

5„ ^ 5(1) = T5„T-i + dtTT-^ (3.2) 

which will, in general, not be a pure g-DO although the i?„ does. However if we suitable 
choose the gauge operator T such that B'^\ as defined by (|3.2|) , is also a purely g-PDO, 
then L^ji) represents a valid new solution to the A^-th g-KdV hierarchy. This is the goal we 
want to achieve in this letter. 



To formulate the DBTs, we follow [|16[ to introduce a g-version of the bilinear potential 
f2(0, t/^) which is constructed from a wave function and an adjoint wave function ip. The 
usefulness of this bilinear potential will be clear when we use it to construct DBTs (see 
below) . 

Proposition 2: For any pair of cf) and if), there exists a bilinear potential n((^, V') 
satisfying 

(a,n(0,V^))=0V^ (3.3) 
dtM<P. ^) = res(9- Vi^nH"') (3-4) 
In fact the bilinear potential fi(0, tp) can be formally represented by a g-integration as 

i](0,^) = (5-VV^). 

Motivated by the DBTs for the ordinary KdV (or Kadomtsev-Petviashvili (KP) 
ITgHTsI) hierarchy, we can construct a qualified gauge operator T as follows 

T^=e{<p,)d,<j)^'=d,-au = (3.5) 

where (pi is a wave function associated with the linear system ( p.24| ). It is not hard to show 
that the transformed Lax operator L^^) is a purely g-PDO with order N and the Lax equation 
( ^.19| ) transforms covariantly, i.e. dt^^L^-^^ = [(L^^^)"^^, L^^^]. The transformed coefficients 
{?4j then can be expressed in terms of {ui} and (pi. On the other hand, for a given generic 
wave function (p ^ (pi (or adjoint function ip) its transformed result can be expressed in 
terms of (pi and itself. 

Proposition 3: Suppose (pi is a wave function of the linear system ( p.24| ). Then the 
gauge operator Ti triggers the following DBT: 
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lW = TiLTr' = a," + iiSLiS,"-' + ■ ■ • + 11?' (3.6) 

= (T,^) = liiMiil, (3.7) 

01 

^« = {{Trri^) = (3.8) 



where L^^), and ?/^(^) satisfy {^J^), and ( p:25| ) respectively. 

Just like the DBTs for the ordinary KdV hierarchy, the DBT triggered by the gauge 
operator Ti is by no means the only transformation in this g-analogue framework. We have 
another construction of DBT by using the adjoint wave function associated with the adjoint 
linear system ( p.25| ). Indeed, for a given adjoint wave function ipi we can construct a gauge 
operator 

s, = ri(v^r^)5,-Vi = id, + A = (3.9) 

which preserves the form of the Lax operator and the Lax equation. 

Proposition 4: Suppose ipi is an adjoint wave function of the adjoint linear system 
( p.25| ). Then the gauge operator 5*1 triggers the following adjoint DBT: 

L(^) = SiLS^^ (3.10) 

0-H..0)^^ (3.11) 

^« = ((5-)*^) = ^^^Mkil, ^^^^p, (3.12) 

^1 



where L'^^\ 0*^^^ and ijj^^'^ satisfy (|2.19|) , (|2.24|) and ( p. 251) , respectively and Wg is obtained 
from Wq by replacing dq with d*. 

So far, we have constructed two elementary DBTs triggered by the gauge operators 
Ti and 5*1. Regarding them as the building blocks, we can construct more complicated 
transformations from the compositions of these elementary DBTs. However, we will see 
that it is convenient to consider a DBT followed by an adjoint DBT and vice versa because 
such combination will frequently appear in more complicated DBTs. So let us compose them 
to form a single operator Ri which we call binary gauge operator. The construction of the 
binary gauge operator Ri can be realized as follows: first we perform a DBT triggered by 
the gauge operator Ti = 9{(l)i)dq(f)i^ and the adjoint wave function tpi is thus transformed 
to ipi'^ = {{Ti^yipi) = —6{(f)i^)9{Q{(j)i,ilJi)). Then a subsequent adjoint DBT triggered by 
Si'^ = 6^'^ {{ip[^^)^'^)d^^ (ipi^^) is performed and the composition of these two transformations 
gives 

R, = S['^T, = 1 - 0l^](0l, (3.13) 

Proposition 5: Let 0i and ■j/'i be wave function and adjoint wave function associated 
with the linear systems ( ^^.241 ) and ( p^.25D , respectively. Then the gauge operator Rx triggers 
the following binary DBT: 
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0(1) = = - 



(3.14) 
(3.15) 
(3.16) 



where L^'^\ 0^^) and satisfy (|]T|), (|]24D and (|]2|) respectively. 

We would like to remark that the construction of the binary gauge operator Ri is in- 
dependent of the order of the gauge operators T and S. If we apply the gauge operator 
5*1 followed by Ti'\ then a direct calculation shows that R = T^^^ Si has the same form as 

The remaining part of this section is to consider the iteration of the DBTs by using the 
DBT, the adjoint DBT, and the binary DBT triggered by the gauge operators T, S, and 
R, respectively. For example, by iterating the DBT triggered by the gauge operator T, we 
can express the solution of the A^-th g-KdV hierarchy through the g-deformed Wronskian 
representation. This construction starts with n wave functions 0i, 02, ■ ■ ■ , 0n of the linear 
system ( p.24| ). Using 0i, say, to perform the first DBT of Proposition 1, then all 0j are 
transformed to 0-^'' = (Ti0j). Obviously, we have 0^^^ = 0. The next step is to perform a 



subsequent DBT triggered by (f)^\ which leads to the new wave functions (f)\'^' with 
Iterating this process such that all the wave functions are exhausted, then an n-step DBT 
triggered by the gauge operator T„ = {dq — a^~^^){dq — a^j^_i^) ■ ■ ■ {dg — ai) is obtained. 



(2) 



d2) 



0. 



{dq(j)^P) / It is easy to see that T„ is an nth-order g-DO of the form 



where a) 

Tn = dg + an-idg~^ -|- ■ ■ ■ -|- ao with Oj defined by the conditions (TnCpj) = 0,j = 1, 2 
Following the Cramer's formula it turns out that Oj = — VF^^*) [0i, ■ ■ ■ , 0„]/l^g[0i, ■ ■ ■ , 0„] 
where W^'^^ is obtained from Wg with its i-th row replaced by (9^ 0i ),■■■, (9^0^). This 
implies that the n-step transformed wave function 0(") (0 7^ 0j) is given by 



(T„ 



Wg[cP, 



(3.17) 



and the n-step gauge operator T„ can be expressed as 



«01 



idq(j)n) dq 



(3.18) 



where it should be realized that in the expansion of the determinant by the elements of the 
last column, dq have to be written to the right of the minors. 

Next let us turn to the iteration of the adjoint DBT. In this case, the n-step gauge 
operator can be constructed in a similar manner by preparing n initial adjoint wave 
functions tpi, ■ ■ ■ ^tpn such that 5"^ = {dq + f3^^-^'>){dq + fit'i^^) ■ ■ ■ {dq + /5i) = d^ + 
YJi=i d^h- Using the required conditions {{S^^^^Yipi) = 0,j = l,---,n we obtain bi = 
■ ■ , ipn]/Wq[ipi, • ■ ■ , 'ipn] which givc the n-step transformed adjoint wave function 



^(") = {{s-'m 



Wqiij 



I5 



Wq[iJi,---,iJn] 



(3.19) 
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and the gauge operator 



(5, 



(3.20) 



Finally, we shall construct n-step binary DBT by preparing n wave functions (pi and n adjoint 
wave functions tpi at the begining. Then we perform the first binary DBT by using the gauge 
operator Ri which is constructed from 0i and tpi as in ( |3.13| ). At the same time, all (pi and ipi 
are transformed to (f)\^^ = (-Ri0i) and ipl^^ = {{Ri^yipi), respectively except (pi'^ = ipi'^ = 0. 
We then use the pair {02^'*, to perform the next binary DBT. Iterating this process 
until all the pairs {(pi, ipi} are exhausted, then we obtain an n-step gauge operator of the 
form Rn = I — Cjdg^ipj- Solving the conditions {Rn(pj) = 0, j = 1, ■ ■ ■ , n, we obtain 
the coefficients Cj = det n^*^ detfi where fij^ = ^l{(pi,-ipj) and fi*^*^ is obtained from n with 



its i-th column replaced by 



This leads to the following representations for R^. 



detn 



^11 ■ ■ ■ ^Ira 01 



(3.21) 



Moreover, the n-step transformed bilinear potential f2(0("\^(")) can be expressed in terms 
of binary-type determinant as 



detfi 



VLxn ^^(01,^) 



^nl ■■■ ^nn il(0„,^) 
fi(0,^l) ■■■ n(0,V'n) ^{CP,^) 



(3.22) 



where 0(") = (i?„0) and ^^"^ = ((i?~^)*V')- 



IV. Q-DEFORMED KDV HIERARCHY {N = 2) 

This section is devoted to illustrate the DBTs for the simplest nontrivial example: q- 
deformed KdV hierarchy (A^ = 2 case). Let 



n = 1, 3, 5, • • • 



L = Uidq + uo 

then the Lax equations 

d,„L=[LT,L], 

define the evolution equations for ui and uq. In particular, for the ti-fiow, we have 

dt^ui = x{q - l)dt^uo 



(4.1) 



(4.2) 



(4.3) 
(4.4) 
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which is nontrivial and recovers the ordinary case as q goes to 1. For higher hierarchy flows 
the evolution equations for Ui and Uq become more complicated due to the non- commutative 
nature of the g-deformed formulation. 

We now perform the DBT of Proposition 3 to the Lax operator (|4.1| ), then the trans- 
formed coefficients become 

u^^^ = e{ui)-ai + e'^{ai) (4.5) 
= 6{uq) + {dgUi) + (g + l)6{dqai) - a^u^ + d{ai)ut'^ (4.6) 

Since 0i is a wave function associated with the Lax operator ( [4.1D , i.e. L0i = Ai0i, one can 
easily verify that {dqai) + 6{ai)ai + Uiai + uq = \i and hence Eqs.(^4.5|) and ( [4.6|) can be 
simplified as 

M«-Mi = x(g-l)(4'^-Mo) (4.7) 
u^Q^ - Uq = dq{ui + ai + e{ai)) (4.8) 

1 /2 

Furthermore, using the facts that dt^(t)i = (L/ 0i) = {dq(f)i) + so(/>i and ui = 9{sq) + sq, we 
can rewrite ( [4.8|) as 

=Mo + 9Aln0i^(0i) (4.9) 

Eqs. (|4.7| ) and (|4.9| ) are just the desired result announced in Sec. III. 

Similarly, by applying the above analysis to the adjoint and binary DBTs, we obtain the 
following result: 

Proposition 6: Let 0i and ipi be (adjoint) wave function associated with the Lax 
operator ( |4.1| ). Then under the DBT, adjoint DBT, and binary DBT, the transformed 
coefficients are given by 



with 

(1) 



u] ' — Ui = x{q — l)('Uo^'' — Mo) (4-10) 



n'=Uo + dgdt. In 01^(01), (DBT) (4.11) 
= uo + dgdt, In^i^-i(^i), (adjoint DBT) (4.12) 
=Uo + dgdt, ln^]n0(^]n), (binary DBT) (4.13) 

Eqs.( 4.11 )-( 4.13|) effectively represent the 1-step transformations. To obtain the n-step DBT, 
adjoint DBT and binary DBT we just need to iterate the corresponding 1-step transforma- 
tions successively by inserting the triggered wave function ( |3.17|) , adjoint wave function 
(|]T|) and bilinear potential (|]22D into the logarithm in Eqs. ( pID , ( CT ) and ( p3l) , 
respectively. 

Proposition 7: Let (pi and ipi {i = 1,2, ■■■,n) be (adjoint) wave functions associated 
with the Lax operator ( |4.1| ). Then under the successive DBT, adjoint DBT and binary DBT 
of Proposition 6, the ra-step transformed coefficients are given by 



in) 



Ui=x{q-l){u^Q^ -uo) (4.14) 
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with 

„W =Uo + dgdt. In ■ ■ ■ , <Pn]0{W,[(l),, ■ ■ ■ , 0„]), (DBT) (4.15) 

= no + dgdt. In lyj^i, ■ ■ ■ , tlJn]e-\Wg[ilJi, ■ ■ ■ , V-n]), (adjoint DBT) (4.16) 

mJ") =uo + dgdt. In det fi^(det fi), (binary DBT) (4.17) 

Eqs.( 4.l5| )-( 4.17| ) provide us a convenient way to construct new solutions from the old ones. 
Especially, starting from the trivial solution {ui = uq = 0) we can obtain nontrivial multi- 
soliton solutions just by putting the (adjoint) wave functions into the formulas ( [4.15D -( ^A7|) . 
For example, the wave functions (pi {i = 1, • • •,n) associated with the trivial Lax operator 
L = d'^ satisfy 

d'^^, = p^cj),, dtj, = (d^cP,) ^ p, (4.18) 

which give the following solutions 

oo oo 

(f)i{x,t) = Eqipix) exp{Y,pf^%k+i) + liEq{-pix) exp(- +H2fc+i) (4.19) 

fc=0 fc=0 

where 7^ are constants and Eg{x) denotes the g-exponential function which satisfies 
dqEg{px) = pEq{px) and has the following representation: 

E,(x)=expEi^^x1 (4.20) 

Substituting ( [4.19| ) into (|4.15|) - ([4.17| ) gives us n-solition solutions. Note that for the soliton 
solutions constructed out from the trivial one {ui = = 0) as described above, they satisfy 
a simple relation: m^"^ = x{q — 1)u"q \ This is just the case considered by Haine and Iliev in 



|Tl|,|T3|. In general, it can be shown |T^ that the solutions of the g-KdV hierarchy can be 
expressed in terms of a single function called tau-function such that 

u,{x,t) = dt, In^^ 



= x{q - l)dqdt, lnr,(x,t)^(r,(x,t)) (4.21) 

Eq. ( ^4.21|) together with Proposition 7 shows that for a given solution (or r^) of the g-KdV 
hierarchy, the transformation properties of can be summarized as follows 

r,^r(") = iyj0i,---,0„]-r, (DBT) 

-r, ^ ^i"^ = d'\Wq[iJi. ■ ■ ■ , ^n]) ■ r, (adjoint DBT) (4.22) 
^ r^(") = det • (binary DBT) 

which implies that the Wronskian-type (or binary-type) tau-function can be viewed as the 
n-step transformed tau-function constructed from the trivial solution (r^ = 1). 
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V. CONCLUDING REMARKS 



We have constructed the elementary DBTs for the q-KdV hierarchy, which preserve the 
form of the Lax operator and are compatible with the Lax equations. Iterated application 
of these elementary DBTs produces new soliton solutions (tau-functions ) of the g-KdV 
hierarchy out of given ones. Following the similar treatment and using the tau-function 
representation u^^i = dt^ hiO^ {jq) /Tq |14] for > 2 we can reach the same result as (^^22]) 



except now dTq/dt-^ = 0. 

In fact these DBTs can be applied to the g-deformed KP hierarchy without difficulty 
by considering the Lax operator of the form L = dq + J2'^oUidq\ The g-KdV is just a 
reduction of the g-KP by imposing the condition {L^)+ = ■ Since the ordinary KP 
hierarchy admits other reductions which are also integrable. Hence it is quite natural to ask 
whether there exist g-analogue of these reductions and what are the DBTs associated with 
them? We hope we can answer this question in the near future. 
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